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fv^ ' The statistical properties of Charmonium energy spectrum determined by the Bethe- 

Salpeter equation are investigated. It is found that the regular motion of the cc system 

> . 
^\ can be expected at a small value of color screening mass but the chaotic motion at a large 

("^ ' one. It is shown that the level mixing due to color screening serves as a new mechanism 

^^ , resulting in the J/ip suppression. Moreover, this kind of suppression can occur before the 
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0^ ' color screening mass reaches its critical value for J/ip dissociation. It implies that a strong 






X 



3/tp suppression is possible in the absence of dissociation of 3/4>- 
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Strongly interacting matter of sufficiently high density is predicted to undergo a transi- 
tion to a state of deconfined quarks and gluons. Deconfinement occurs when color screening 
shields a given quark from the binding potential of any other quarks or antiquarks. Bound 
states of heavy quarks, such as the J/^/" or the T, whose radii are much smaller than those of 
the usual mesons and nucleons, can survive in a deconfined medium until the temperature 
or density becomes so high that color screening prevents their tighter binding. It was theo- 
retically predicted that a suppression of J/i/) or tp production due to the dissociation can be 
found in relativistic heavy ion collision, which can serve as a clear signature for the formation 
of Quark-Gluon-Plasma[l]. Subsequently this suppression was observed by NA38 coUabora- 
tion[2]. However, successive research pointed out that such suppression could also exist in 
hadronic matter, even though caused by a completely different mechanism [3]. Recently, an 
anomalously strong J/tp suppression has been observed by the NA50 collaboration[4] and 
there have been a number of attempts(see for example Refs.[5-7]) to explain it. However, 
the mechanism of J/i/> suppression is still a question of debate. 

It is known that, for a complex system, the statistical properties of the quantum spectrum 
and the complexity of the eigenfunctions reveal frequently new features of the system. Be- 
cause of the complexity of the interaction and the temperature and density effect, the heavy 
qq systems in hot matter can be regarded as complex systems. Nevertheless, the spectral 
statistical properties and the complexity of eigenfunctions of heavy qq systems have not yet 
been discussed. 

In this paper we study the spectral statistical properties and the complexity of eigenfunc- 
tions of the cc system. The influence of color screening on the statistical properties and 
the mixing of eigenfunctions which characterizes the complexity of eigenfunctions will be 
examined. And an alternative interpretation of the strong J/t/j suppression is proposed. 

To investigate the statistical properties of the cc system a reasonable energy spectrum 
is crucial. Usually the bound states of the cc system were studied in the nonrelativistic 
formulation[8,9]. However, the motion of a quark and an antiquark in a meson should be 
relativistic. As pointed out in Ref. [10], the kinetic energy of the cc system is about 13% of 
the total energy and the ratio of relativistic corrections to the quark mass will not decrease 
with the increase of quark mass. As a result the bound state equation for J/tp should be 
relativistic. Since the Bethe-Salpeter (BS) equation[ll] is the unique effective relativistic 



equation of two-body bound states and is consistent with quantum field theory, we take the 
BS equation to determine the energy spectrum of charmonium in this paper. 

It is known that the full bound state BS equation, written in the two-sided notation, 
reads [12] 



where rji = ^ ""'^ (i=l,2), xp is the wave function for the quark-antiquark system with total 
four momentum P, p is the relative four momentum. V is the interaction kernel that acts on 
Xp and formal products VXpip') take the form V Xp{p') = K Xp(p') + 7m ® 7'^K Xp{p')j 
in which Vs and Vv are the scalar, vector potential respectively. 

With the standard reduction and spin-independent treatment being implemented, the 
three-dimensional spin- independent reduced Salpeter equation can be written as[13], 

(M - £i - S2)0(p) = / -^ V Frip,p')v.i\p- p'lmP) (2) 



where (f){p) = Jdp'^ Xp{p'\p) is the three-dimensional equal-time BS wave function. M is 
the mass of the qq bound state. Ei = (p^ +m.^)2 ^ i=l,2 represent a quark and an antiquark 
respectively. The functions _F^* and F" are the transformation coefficients of the spin- 
independent vector, scalar potential respectively when the four dimensional BS equation 
is reduced to the three dimensional spin-independent BS equation. Their expressions have 
been given explicitly in Ref. [13]. 

It is definite that the momentum dependence of the interaction is treated exactly in Eq.(2). 
And Eq.(2) gives a well-defined eigenvalue problem for the masses M of the qq bound states 
in the momentum space. To solve Eq.(2), we take the color screening potential between the 
quark and antiquark[8]. In momentum space the color screening potential is written as 

In our calculation, the reduced BS equation which has been described well in Refs. [12,13] 
and a color screened quark-antiquark potential[13] are implemented. In a momentum space, 
the potential with the color screening can be written as 

Vs{\p-p'\) = y/{p-p') - ^ I , (3) 

and 

K(|P-P'|) = -^-— ^— -, (4) 

37r^ [(p - p') + M ] 



where Vs, Vv refer to the scalar, vector potential respectively, and g is the string tension, Qs 
the effective coupling constant, and ^ the Debye screening mass. The parameters including 
the value of the /i at zero temperature (denoted as /io hereafter) are determined by a least 
square fitting to the experimental data of IS, 2S and 3S states of the cc system. The 
fixed parameters are a=0.22 GeV^, /io=0.06 GeV, 7nc=1.474 GeV, as{cc) — 0.47, which are 
within the scope of customary usage. Some of the obtained spin-averaged energy (mass) 
[13] spectrum of the cc bound states at zero temperature (T — 0) are listed in Table I. From 
Table I, one can observe that the mass spectrum obtained at present is more consistent with 
the experiment than the previous ones[8,9]. With the parameters a, rric and as{cc) fixed 
above, we evaluate the mass spectrum of the cc bound states for different ^, so that the 
temperature dependence of the spectrum can be investigated. The calculation indicates that 
the critical value of the screening mass for J/i/) to be dissociated is about ^c=0.90 GeV. 

As the spectrum has been determined, we study the spectral statistics. To this end we ex- 
ploit the well established random matrix theory (RMT)[14,15]. The RMT was initially used 
to describe the statistical properties of complex nuclear spectra in the early of 1950's[14]. 
In the RMT, the matrix of the Hamiltonian of the system, which is sufficiently complex and 
(or) governed by a very complicated or even unknown dynamics, is replaced by a random 
matrix with well-defined symmetry properties. Any desired quantity of the system is de- 
termined by performing an average over an ensemble of the random matrices of which the 
elements are distributed according to a certain probability distribution. The distribution is 
usually taken to be a Gaussian with a parameter d, 

P{H) oc exp(-Tr{HH+}/2cP), (5) 

and the corresponding random matrix ensemble is called Gaussian ensemble. The under- 
lying space-time symmetries obeyed by the system impose some important restrictions on 
the admissible ensemble. If the Hamiltonian is time-reversal and rotational invariant, the 
Hamiltonian matrices can be chosen to be real symmetric. The corresponding ensemble is 
called the Gaussian orthogonal ensemble(GOE). If the Hamiltonian is not time-reversal in- 
variant, irrespective of its behaviour under rotations, the Hamiltonian matrices are complex 
Hermitian. The corresponding ensemble is called the Gaussian unitary ensemble (GUE). 
If the system is time-reversal invariant but not invariant under rotations, and it has half- 
odd-integer total angular momentum, the matrices are quaternion real. The corresponding 



ensemble is called the Gaussian symplectic ensemble (GSE). 

If the Hamiltonian possesses other invariances, the corresponding ensemble can be the 
Gaussian unitary ensemble (GUE), or the Gaussian symplectic ensemble (GSE). The RMT 
has been applied successfully in many areas of physics, such as nuclear physics [16], quantum 
field theory [17] and the physics of disordered systems[18]. The success of the RMT is rooted 
in the fact that there exist so-called universal quantities which do not depend on the details 
of the dynamics but only on the underlying symmetries. This enables one to separate generic 
features from properties that do depend on the details of dynamics. The universal quantities 
usually mean the nearest neighbor level spacing distribution P{s) and the spectral rigidity 
A3 (I/) which are powerful to analyze the spectral statistical properties of complex systems. 

In present case the cc system is time-revsersal and rotational invariant. It is then a member 
of the GOE and its spectrum should be analyzed by comparing with the predictions of the 
GOE. To make comparison with the universal and dimenssionless results of the RMT, it 
is necessary to make a transformation or normalization on the spectrum. This operation 
is called the unfolding. The procedure is as follows. For a spectrum {£"*}, one separates 
its smoothed average part from its fluctuating part at first, then counts the number of the 
levels below E so that one can define a staircase function N{E) of the spectrum, 

N{E) = Na4E) + Nfi^rAE). 
The unfolded spectrum can be finally obtained with the mapping 

Ei=Nav{E,). 

This unfolded levels Ei are obviously dimensionless and have a constant average spacing of 
one, but the actual spacings exhibit frequently strong fiuctuations. The nearest neighbor 
level spacings are defined as Si — Ei+i — Ei. The distribution P{s) is defined as that P{s) ds 
is the probability that s; lies within the infinitesimal interval [s,s + ds]. 

It has been shown that the nearest neighbor level spacing distribution P{s) measures the 
level repulsion (the tendency of levels to avoid clustering) and short-range correlations be- 
tween levels, which mean that, for a given level, its repulsion to the neighbor levels is limited 

within several level spacings[18,19] . For the GOE, the distribution can be expressed explicitly 

2 
as P{s) = f- s exp( — ^^-^), which can be deduced from Eq.(5) [15]. With the Brody param- 
eter w in the Brody distribution P„ (s) = (l-|-a))as'^exp(-QS^+")(a = [r(fi^)]^+",r(a-) is 



the F- function), the transition from regularity to chaos can be measured quantitatively. It 
is evident that cu — 1 corresponds to the GOE distribution, while a; = to the Poisson-type 
distribution. A value < a; < 1 means an interplay between the regular and the chaotic. 
As to the spectral rigidity A3(L), it is defined as 

1 r^/^ 

As{L) = (min^ ^ - / [N{x) - Ax - Bfdx), (6) 

^ J-L/2 

where A''(a;) is the staircase function of a unfolded spectrum in the interval [—L/2,x]. The 
minimum is taken with respect to the parameters A and B. The average denoted by {• • •) 
is taken over a suitable energy interval over x. Thus from this definition A3(L) is the local 
average least square deviation of the staircase function N{x) from the best fitting straight 
line. For the GOE the expected value of A3(L) can only be evaluated numerically, but it 
approaches the value A3(L) = -\ {InL — 0.0687) for large L. It has also been shown that 
the spectral rigidity A3(I/) signifies the long-range correlations of quantum spcctra[18,19] 
which make it possible that for a chaotic spectrum very small fluctuation of the staircase 
function around its average can be found in an interval of given length (the interval may 
cover dozens of level spacings). 

It has been generally accepted that if the spectral statistical properties of a quantum 
system approach to the GOE, the motion of the system is chaotic, and if the statistical 
properties appear as Poisson-types {P{s) — exp( — s) and A3(L) = j^) ), the motion is 
regular [18]. 

To perform a meaningful RMT analysis one has to sort the spectrum in symmetry sectors 
corresponding to the symmetry of the Hamiltonian. Since the symmetry invariant subspaces 
are orthogonal to one another, each such symmetry invariant subspace is characterized by a 
specific set of quantum numbers. The RMT analysis is performed on sets of the eigenlevels 
having the same quantum numbers. In view of the above consideration we take the first 
200 levels out of 300 s-wave (spin-weighted) eigenvalues which have the same quantum 
numbers (e.g., / = 0). We have noticed that p-wave eigenvalues (spin- weighted) have the 
same statistical properties as those of the s-wave ones. The spacing distribution P{s) and 
spectral rigidity A3 for several color screening masses ji are shown in Fig.l. The figure 
shows obviously that, when /i is small, for instance 0.06 GeV, both the P(s) and A3 display 
the Poisson-type distribution. This indicates that the spectrum is regular and the levels 
are uncorrelated. The motion of the cc system is then regular. With the increase of ji (for 



example /i = 0.20 GeV), the P{s) and A3 depart from the Poisson- types. Then the spectrum 
becomes less regular and the level correlation and repulsion increase gradually. As /i is large 
enough, say /i=0.60 GeV, both the P{s) and the A3 appear the GOE behavior. This means 
that the spectrum becomes extreme irregular (chaotic) and the levels are strongly correlated. 
Strong level repulsion can thus be expected and the motion of the cc system is chaotic. 

Fitting to the Brody distribution, we get the Brody parameter co = 0.22,0.51,0.89 for 
/i = 0.06, 0.20, 0.60 (GeV) respectively. This also indicates that a transition from regular 
to chaotic motion can be induced in the cc system by the increasing of the color screening 
and the chaotic motion takes place when the hot matter is of high temperature. Recalling 
the idea that the P{s) and A3 measures the correlation among the levels[18,19], one can 
recognize that the correlation of the cc resonances can be manifested by the P{s) and A3. 
The calculated result shows thus that the cc resonances are strongly correlated at a high 
temperature. This information may be useful for us to analyze the relevant experiments. 

As a matter of fact, both the level correlation and level repulsion stem from level mixing 
[20] which characterizes the complexity of eigenfunctions. We analyze then the influence of 
the color screening on the level mixing. For this purpose, we take the s-wave eigenfunctions 
of the cc system at temperature zero |A;) (J/i/", tp , i>" , ■ ■ ■, /i=0.06 GeV, equivalent to T=0) 
as a basis and expand the eigenfunctions at a finite temperature \i) in this basis as 



|i)=^cL|fc), (7) 



k 



With Efe(4)' = 1- 

In order to measure the level mixing and delocalization of eigenfunctions at a finite tem- 
perature, we define the following probability 

where c^ is the probability for the eigenfunctions at zero temperature to appear. AN{c^) 
is the number having the probability between c^ and c^ + 5 with 5 being a small interval 
of c^. Nt is the total number of the eigenstates involved (at present Nt — 200). For the 
J/ip state (the ground state) at a finite temperature, the P{c^) versus c^ for a few values 
of n are shown in Fig. 2. It is evident that, when fi is small, for example 0.10 GeV, the 
components with nonzero probabilities are limited within the lowest states. It means that 
the components with small probabilities prevail and the distribution width of Pic^) is rather 



small. In other words, the level mixing is localized. With the increase of /i, the number 
of the components with nonzero probabilities goes up. This leads to the increase of P{c^) 
in the region of large probabilities (see the case ^=0.20 GeV). As /i approachs to fic, the 
components with nonzero probability emerge on a large scale and become rather random (as 
shown by the cases /i=0.60, 0.90 GeV). Then, a wide distribution of P(c^) appears and the 
level mixing is strong and delocalized. Referring to Fig.l, one can find that the quantum 
chaotic motion occurs when a strong level mixing emerges. It can also be found that, for 
higher eigenstates at a finite temperature, the level mixing changes with ^ in a similar way 
as J/V> state. Furthermore the strong level mixing appears for these eigenstates even if /i 
takes a smaller value. 

The feature of the effect of the level mixing induced by the variance of jj, is that the 
probability for the J/i/) state at zero temperature to appear in the J/tp state at a finite 
temperature descends with the increasing of /i. Present result shows that at /i=0.10, 0.20, 
the probability is about 90%, 70%, respectively. As /i increases, the probability decreases 
distinctly. For instance, when jj,— 0.60 GeV, the probability goes down to 10%. Referring 
to Eq.(7) one can know that the neighboring levels (e.g. tp , tp ...) at a finite temperature 
will also contribute to the J/i/> at zero temperature. However, their production probability 
is much smaller than that of J/?/> at a finite temperature (e.g. in p-p collisions, the ratio 
of V' to J/ip is about 1/7[21], and in Pb-Pb collision it is predicted to be 4%[22]). Then 
their contribution to the J/?/> at zero temperature should be very small. One can hence 
easily recognize that the production of J/xp at zero temperature, which means in fact the 
probability for the charmonium states at a finite temperature to include the J/i/> state at 
zero temperature, will be suppressed with the growth of /i and a strong suppression can 
be expected when /i is large enough. On the other hand. It has been well known that, in 
the deconfinement scheme, the suppression due to the dissociation of 3/ip occurs only if /i 
reaches or exceeds the critical value Hc [1]. Present analysis indicates that, the suppression 
due to the level mixing can take place before /i reaches the critical value for the dissociation 
to happen. In addition, we find that for p-wave energy spectrum of the cc system, the 
transition to the chaos and the suppression due to the level mixing occur too. Moreover, 
the similar behavior can be found in the bb system. 

To apply the above discussion to an actual physical situation, one needs to know the spe- 



cific dependence of the Debye screening mass /i on the temperature T. Even though several 
studies on the properties of QCD at finite temperatures and densities have been accomplished 
(see for example Refs. [23-27]) and the general feature of the temperature dependence fi{T) 
has been well established, the explicit expression fi{T) has not yet been determined uniquely 
because a temperature-dependent running coupling constant is involved. Then the density 
dependence of /i is not clear now. Any way, it is certain that the /i(r) increases with T 
increasing. The above calculated result manifests that, as the statistical properties of the 
mass spectrum of the cc system is taken into account, the critical value Hc for the strong 
suppression to take place gets lowered obviously from that in the dcconfinement scheme 
(from 0.9 GeV to 0.6 GeV). Then the corresponding critical temperature Tc decreases too. 
Consequently, much stronger suppression than that predicted in the deconfinement scheme 
can emerge in the same experimental status. The level mixing, or the chaotic motion, due to 
the color screening is then an alternative mechanism to induce the strong J/xp suppression. 
In conclusion, in this letter we have investigated the level statistics of the Charmonium 
energy spectrum in the framework of the BS equation and found that, with the increasing 
of the color screening mass, the transition from regular to chaotic motion can be induced for 
the cc system. The level mixing due to the color screening is then an alternative mechanism 
resulting in the strong J/ip suppression, and furthermore this kind of suppression can occur 
before /i reaches its critical value for the dissociation of cc. This implies that a strong J/ip 
suppression is possible in the absence of dissociation of the J/V'- 
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Table I. The calculated spin-averaged mass spectra M„; of the cc bound states and the comparison with experimental 
data and previous calculations. 
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Figure Captions 



Fig.l The spacing distribution P{s) and the spectral rigidity A3 of the s-wave Charmonium energy spectrum at 
different values of color screening masses jj, measured in GeV. In (a) the dashed-line for Poisson distribution, 
the solid one for the GOE distribution, and the histograms for our numerical results. 

Fig. 2 The probability density P(c^) of -J/tp versus c^ for various values of fi. The solid histogram for /i=0.90 GeV, 
the dotted one for fj,=0.60 GeV, the dash-dotted one for jj.=0.20 GeV and the dashed one for /i=0.10 GeV. 
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